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Abstract 

An  opposed  jet  configuration  featuring  fractal  generated  turbulence  was  used  to  characterise  premixed 
fuel  lean  JP-10  flames  with  equivalence  ratios  (0)  from  0.20  to  0.80  stabilised  against  combustion  products 
with  a  temperature  of  1720  K.  The  cold  flow  turbulent  Reynolds  (Re,)  number  was  set  to  122  and  the  esti¬ 
mated  Damkohler  (Da)  numbers  varied  from  below  0.3  (0  =  0.20)  to  6  (0  =  0.80)  with  laminar  flame 
parameters  obtained  from  computations  of  the  corresponding  opposed  jet  configuration.  The  gradual  tran¬ 
sition  from  a  turbulent  flame  to  a  supported  burning  mode  related  to  the  Homogeneous  Charge  Diffusion 
Ignition  (HCDI)  flameless  oxidation  regime  was  detected  with  features  typical  of  conventional  flames  grad¬ 
ually  disappearing  for  0  <  0.60.  A  second  mode  transition  to  the  corrugated  flamelet  regime  at  0  =  0.80  is 
shown  to  coincide  with  a  gradual  change  from  gradient  to  counter-gradient  turbulent  transport.  A  density 
segregation  based  technique  was  applied  to  Particle  Image  Velocimetry  (PIV)  images  to  detect  the  instan¬ 
taneous  location  of  reaction  zones  and  to  provide  conditional  and  unconditional  statistics  to  illustrate  the 
combustion  regime  transitions.  Flame  surface  area  characteristics  and  turbulent  burning  velocities  were 
also  derived  to  further  characterise  the  burning  modes. 

©  2012  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 
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1.  Introduction 

Weinberg  [1]  and  Hardesty  and  Weinberg  [2] 
suggested  the  introduction  of  heat  exchangers  that 
use  combustion  products  to  preheat  reactants  and 
introduced  the  term  Excess  Enthalpy  Combustion 
(EEC).  The  result  is  an  increase  in  the  overall 
efficiency  and  combustion  chamber  temperatures 
could  be  decreased  without  flame  extinction 
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leading  to  a  reduction  in  the  formation  of  oxides 
of  nitrogen.  High  temperature  air  combustion 
(HiTAC),  using  the  concept  of  exhaust  gas  recir¬ 
culation  (EGR),  was  subsequently  characterised 
by  Katsuki  and  Hasegawa  [3]  and  the  term  flame¬ 
less  oxidation  introduced  by  Plessing  et  al.  [4]  and 
Wiinning  and  Wiinning  [5],  Cavaliere  and  de 
Joannon  [6]  summarised  progress  in  combustion 
technologies  in  terms  of  reaction  stabilisation  via 
preheating.  Dally  et  al.  [7,8]  investigated  the  char¬ 
acteristics  of  Moderate  and  Intense  Low  Oxygen 
Dilution  (MILD)  combustion  in  a  laboratory 
scale  furnace  as  well  as  in  a  hot  coflow  burner. 
Mastorakos  et  al.  [9]  used  the  counterflow  geome¬ 
try  to  stabilise  ultra  lean  premixed  methane  flames 
against  hot  combustion  products  in  the  opposing 
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stream.  It  was  found  that  when  the  temperature  of 
the  hot  product  stream  exceeded  1 550  K,  extinc¬ 
tion  was  not  observed  even  at  an  equivalence  ratio 
of  0.20.  Furthermore,  the  residual  oxygen  mole 
fraction  in  the  hot  products  was  found  to  have 
no  significant  impact  on  flame  stabilisation.  Gei- 
pel  et  al.  [10]  used  a  counterflow  configuration 

[11]  to  further  characterise  the  stability  of  lean 
premixed  methane  flames  burning  against  hot 
combustion  products  at  temperatures  between 
1520  K  and  1820  K.  The  related  Homogeneous 
Charge  Diffusion  Ignition  (HCDI)  regime  was 
identified  computationally  by  de  Joannon  et  al. 

[12]  using  a  laminar  opposed  jet  configuration. 
The  HCDI  regime  is  of  relevance  to  lean  premixed 
pre-vapourized  (LPP)  Gas  Turbine  applications 
[12]  and  of  particular  interest  under  ultra-lean 
conditions  where  flame  stability  issues  come  to 
the  fore  [13],  The  current  work  quantifies  the 
combustion  mode  transition  from  a  propagating 
premixed  turbulent  flame  to  a  HCDI  related 
flameless  oxidation  regime  through  the  applica¬ 
tion  of  Particle  Image  Yelocimetry  (PIV),  a  den¬ 
sity  segregation  technique  [14]  and  the  OH-LIF 
technique  of  Kerl  et  al.  [15],  The  study  expands 
upon  past  work  by  providing  (i)  flow  field  statis¬ 
tics,  (ii)  conditional  statistics,  (iii)  scalar  fluxes 
and  (iv)  novel  data  on  flames  burning  JP-10.  (v) 
The  influence  of  turbulence  on  the  combustion 
regimes  is  analysed  using  a  Damkohler  number 
(Da)  based  approach  [16], 


2.  Experimental  techniques 

The  burner  layout  was  similar  to  that  used  by 
Geyer  et  al.  [17]  and  Bohm  et  al.  [18]  with  two  sig¬ 
nificant  modifications  as  shown  in  Fig.  1.  Firstly, 
the  lower  nozzle  was  fitted  with  a  perforated  plate 
(hole  size  =  4  mm;  blockage  ratio  (BR)  =  45%) 
used  to  stabilise  lean  (cj>  =  0.75)  premixed  CH4 
flames  50  mm  upstream  of  the  nozzle  exit  in  order 
to  generate  hot  combustion  products.  Secondly, 
the  upper  nozzle  was  fitted  with  the  fractal  grid 
[19-21]  previously  used  by  Geipel  et  al.  [11]  in  a 
parametric  study  of  turbulence  generation  in  iso¬ 
thermal  flows.  The  optimal  fractal  grid  [11]  had  a 
BR  of  65%  with  maximum  and  minimum  bar 
widths  of  2  mm  and  0.5  mm  as  shown  in  Fig.  2 
along  with  the  conventional  grid  used  in  the  lower 
nozzle.  The  nozzle  separation  was  set  to  one  nozzle 
diameter  D  (=30  mm).  The  flames  stabilised  on  the 
grid  in  the  lower  nozzle  featured  a  bulk  velocity 
(14)  of  1  m/s  (at  298  K)  resulting  in  ~5%  oxygen 
by  volume  remaining  in  the  products  according 
to  laminar  flame  calculations  [22].  Temperatures 
at  the  lower  nozzle  exit  ( Th)  were  measured  at 
1720  K  using  a  type  R  thermocouple  connected 
to  an  Omega  HH806AU  unit  indicating  heat  losses 
around  6%  to  the  burner  assembly.  Mixtures  of  JP- 
10  (b.p.  460  K  at  1  atm  [23])  and  air  were  intro- 


Fig.  1.  Experimental  configuration.  JP-10/air  mixtures 
introduced  in  the  upper  nozzle  and  hot  products  from  a 
CH4/air  flame  ( </>  =  0.75)  in  the  lower  nozzle. 


duced  in  the  upper  nozzle  at  cf>  —  0.20,  0.40,  0.60 
and  0.80  at  Ub  =  3.5  m/s  (at  298  K)  with  dry  air 
used  as  a  reference.  The  integral  length  scale  of  tur¬ 
bulence  (L/=  3.1  ±  0.1  mm)  was  determined  by 
Geipel  et  al.  [11],  The  JP-10  mixtures  were  heated 
to  473  K  resulting  in  nozzle  exit  temperatures 
Ta  ~  400  K  which  prevented  recondensation.  The 
flow  control  system  was  similar  to  that  used  previ¬ 
ously  [11]  with  an  uncertainty  in  the  flow  rate 
<0.8%  for  each  fluid.  Flow  rates  of  JP-10  were 
metered  and  vapourised  using  Bronkhorst  Cori- 
Flow  M53  and  CEM  W-303A  units.  GC-MS  mea¬ 
surements  performed  on  the  reactant  stream  using 
an  Agilent  7890A  series  GC  with  a  5975C  inert 
MSD/DS  Turbo  El  Bundle  equipped  with  60  m 
DB-1  column  confirmed  that  no  cracking  of  the 
fuel  had  taken  place.  Coflow  velocities  were  set 
to  1  m/s  in  the  lower  nozzle  to  prevent  overheating 
and  0.3  m/s  in  the  upper  nozzle  to  remove  any 
large  scale  bulk  motion  [11].  Velocity  statistics 
was  measured  with  PIV  using  a  120  mj  Solo-New 
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Fig.  2.  Left:  Grid  used  to  stabilise  laminar  flames  in  the 
lower  nozzle.  Right:  Fractal  grid  with  largest  and 
smallest  fractal  bar  widths  of  2.0  mm  and  0.5  mm  used 
in  the  upper  nozzle. 


Wave  Nd:YAG  laser  with  the  upper  stream  seeded 
with  aluminium  oxide  particles  of  about  3  pm 
diameter  [11,24]  with  images  processed  using 
decreasing  interrogation  window  sizes  from 
128  x  128  via  64  x  64  to  32  x  32  with  a  50%  over¬ 
lap  resulting  in  velocity  vectors  spaced  about 
0.4  mm  apart.  Laser  pulses  were  separated  by 
40  ps  to  minimise  spurious  vectors. 


3.  Postprocessing  techniques 

Velocity  statistics  were  obtained  from  1000 
images  for  each  case  using  a  purpose  written 
FORTRAN  program  [24],  A  density  segregation 
algorithm  [14]  was  used  to  detect  isocontours 
marking  the  instantaneous  reaction  zone  (inter¬ 
face)  location.  Such  methods  have  been  used  in 
the  past  [14,24-26]  and  found  to  be  accurate  for 
the  detection  of  conventional  flame  surfaces  with 
statistical  information  readily  generated.  An  anal¬ 
ysis  has  been  presented  by  Steinberg  et  al.  [27]. 
The  accuracy  of  the  method  in  the  presence  of 
broadening  of  the  local  interface  between  hot 
and  cold  streams  was  initially  evaluated  by  con¬ 
sidering  images  of  the  spatial  distribution  of  the 
mean  intensity  (Iimg)  from  all  sets  of  PIV  images. 
The  mean  intensity  in  each  image  was  normalised 
using  the  mean  intensity  in  a  reference  window 
Iwi„  in  the  same  image  where  only  reactants  exist. 
Particle  images  have  significant  deviations  in  Mie 
scattering  intensities  and  it  was  not  possible  to 
obtain  second  order  statistics  directly.  The  ability 
of  the  algorithm  to  detect  isocontours  was  first 
assessed  by  determining  the  probability  of  detect¬ 
ing  hot  products  from  the  lower  nozzle  for  the 
limiting  case  of  no  added  heat  release  (air  emerg¬ 
ing  from  the  upper  nozzle).  The  normalised  inten¬ 
sity  signal  (1  —  Iimg /I Wi„)  represents  the  spatial 
distribution  of  mean  densities  relative  to  that  of 
the  unreacted  mixtures  in  the  upper  stream.  Sub¬ 
sequently,  the  binary  images  detected  by  density 
segregation  were  binned  into  windows  of  16  pixels 
(~0.4  mm)  wide.  The  mean  values  of  these  images 
were  then  calculated  for  each  set  to  derive  the  spa¬ 
tial  distribution  of  the  mean  progress  variable  (c) 


-  directly  related  to  the  overall  mixing  layer  (tur¬ 
bulent  flame  brush)  thickness  (<5f).  Results 
obtained  using  both  techniques  are  shown  for  all 
five  cases  in  Fig.  3.  No  significant  discrepancies 
between  the  mean  spatial  distribution  of  density 
and  the  progress  variable  derived  from  density 
segregation  are  observed. 

The  PIV  images  can  also  be  used  to  estimate  the 
evolution  of  the  interface  thickness  ( d'f )  with  a 
reduction  in  the  Damkohler  (Da)  nurriner  [16], 
The  centreline  profiles  of  the  Mie  scattered  images 
were  shifted  such  that  the  detected  iso-contours 
were  located  at  a  datum  xshift  =  0.  The  shifted  pro¬ 
files  were  normalised  using  the  corresponding 
mean  image  intensity  in  the  reactant  stream  from 
the  upper  nozzle.  Profiles  with  more  than  one  reac¬ 
tion  zone  were  rejected  at  a  rate  between  4.5%  and 
11%  depending  on  the  case.  The  mean  reaction 
zone  density  can  be  mapped  relative  to  the 
detected  isocontours  and  dj  was  obtained  by  cal¬ 
culating  the  mean  of  the  displaced  and  normalised 
signal  (I„orm)  for  every  set  of  1000  images.  Results 
are  shown  in  Fig.  4  and  discussed  in  Section  4. 

The  current  technique  was  further  evaluated  by 
introducing  dry  air  and  methane-air  mixtures 


Fig.  3.  The  impact  of  interface  broadening  on  c 
obtained  using  density  segregation.  Lines:  The  mean 
image  intensity  normalised  by  the  mean  image  intensity 
in  the  reference  window  (1  —  7ims/7wi„).  (O)  The  resolved 
c  when  aligned  with  velocity  vector  locations.  Top  to 
bottom:  Isothermal  case,  <j>  =  0.20,  0.40,  0.60  and  0.80. 


3314 


K.H.H.  Goh  et  al.  I  Proceedings  of  the  Combustion  Institute  34  (2013)  3311-3318 


Fig.  4.  The  evolution  of  the  interface  layer 
thickness.  Solid  line:  Normalised  signals  7„ora,  centred 
on  detected  reaction  zones.  Dashed  line:  Error  function 
fits.  Dot-dashed  line:  c  =  (T  —  T0)/{Th  —  T„)  from  a 
laminar  flame  solution.  Dotted  line:  The  corresponding 
heat  release.  Cases  as  for  Fig.  3. 


with  (j>  =  0.40,  0.60  and  0.80  at  the  upper  nozzle 
and  performing  simultaneous  PIV  and  OH-LIF 
measurements  using  a  Ba(N03)2  (barium  nitrate) 
crystal  based  Raman  laser  configuration  [15]  to 
produce  laser  light  at  a  wavelength  of  281.7  nm. 
The  mean  of  the  OH-LIF  images  was  resolved 
for  each  set  and  subsequently  normalised  with 
the  mean  signals  from  the  exit  streams  of  each 
nozzle  to  obtain  a  reaction  progress  variable 
C  =  (7oh  -~IoH,r)/{IoH,p  -^on,r),  where  )-  and  p 
represent  the  streams  from  upper  and  lower 
nozzles.  The  normalised  signals  from  PIV  and 
OH-LIF  were  compared  and  an  example  is  shown 
in  Fig.  5.  There  are  no  notable  differences  between 
the  different  methods  for  ^  ^  0.6.  The  OH-LIF 
derived  profile  for  4>  =  0.8  includes  a  peak  indicat¬ 
ing  an  elevated  combustion  temperature  consis¬ 
tent  with  the  transition  to  a  conventional 
propagating  premixed  turbulent  flame. 


4.  Results  and  discussion 

Velocity  statistics  along  the  burner  radius 
obtained  2  mm  from  the  upper  nozzle  exit  are  pre¬ 


Fig.  5.  Application  of  the  density  segregation  technique 
and  simultaneous  PIV  and  OH-LIF.  (□)  Mean  PIV 
image  intensity  normalised  by  mean  image  intensity  in  a 
reference  window  (1  —  7, mg //„•„).  (O)  Mean  intensity  of 
normalised  PIV  images  (1  -  {Iimg,iflmn,i)),  where  i  rep¬ 
resents  the  instantaneous  image.  (O)  Detected  mean 
progress  variable.  (A)  Normalised  OH-LIF  signals 
((7 OH  -IoH,r)/(IoHj,-IoH,r)),  where  r  and  p  represent 
reactant  streams  from  the  upper  and  lower  nozzles 
respectively.  Top  to  bottom:  Isothermal  case,  <j>  =  0.40, 
0.60  and  0.80. 


sented  in  Fig.  6  with  the  <j>  =  0.20  and  0.40  data 
sets  omitted  due  to  negligible  deviations  from 
the  isothermal  case.  However,  significant  effects 
on  mean  radial  velocities  were  noted  for  the  </> 
=  0.60  and  0.80  cases.  Velocity  statistics  obtained 
using  conventional  PIV  along  the  stagnation  point 
streamline  can  be  found  in  Fig.  7,  where  the  mean 
axial  velocities  (I/)  and  axial  («')  and  radial  ( v ') 
rms  fluctuations  are  shown.  The  data  suggests  a 
transitional  behaviour  with  no  significant  devia¬ 
tions  in  velocity  statistics  from  the  isothermal  case 
until  <j>  >  0.4  and  with  a  clear  difference  emerging 
at  <j)  =  0.8.  A  significant  amount  of  statistics  was 
obtained  using  the  density  segregation  method. 
As  outlined  above,  the  method  is  well  established 
for  thin  reaction  zones  and  can  be  expected  to 
provide  reasonably  accurate  data  in  the  presence 
of  broadening  of  the  interface  layer.  The  changing 
nature  of  the  combustion  mode  is  apparent  from 
the  progress  variable  statistics  (c  and  c'c')  shown 
in  Fig.  8.  The  mean  location  of  the  detected  reac¬ 
tion  zones  in  the  upper  stream  show  no  significant 
deviations  from  the  isothermal  case  until  tj>  >  0.40. 
The  result  confirms  the  absence  of  a  propagating 
turbulent  flame.  Sample  PIV  images  are  shown 
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Fig.  6.  Normalised  mean  axial  (top  left)  and  radial  (top 
right)  velocities  along  with  axial  (bottom  left)  and  radial 
(bottom  right)  fluctuations  measured  2  mm  from  the 
upper  nozzle  exit.  Isothermal  ( x),  </>  =  0.60  (A)  and  0.80 
(□)  cases  are  shown.  Data  for  </>  =  0.20  and  0.40  have 
been  omitted  as  there  is  no  significant  deviation  from  the 
isothermal  case. 


Centreline  [x/DJ  Centreline  [x/D] 


Fig.  7.  The  mean  axial  velocity  (O)  (left)  and  the  axial 
(□)  and  radial  (A)  velocity  fluctuations  (right)  along  the 
burner  centre  axis.  Top  to  bottom:  Isothermal  case,  tj> 
=  0.20,  0.40,  0.60  and  0.80. 


in  Fig.  9.  Pockets  of  JP-10  combustion  products 
are  starting  to  be  observed  at  4>  —  0.4  and  indicate 
the  gradual  transition  to  conventional  premixed 
flames  as  <f>  — >  0.8.  The  temperature  of  the  hot 
products  from  the  lower  nozzle  is  well  above  that 
of  auto-ignition  and  the  resulting  burning  mode 
for  leaner  mixtures  is  related  to  the  HCDI  regime 
[12]  with  fuel  conversion  increasingly  sustained  by 
the  heat  from  the  combustion  products  [9,10],  The 


Fig.  8.  First  (top)  and  second  (bottom)  moments  of  the 
reaction  progress  variable  (c)  along  the  centre  axis  of  the 
burner  for  <£  =  0.20  (O),  0.40  (O),  0.60  (A)  and  0.80  (□), 
along  with  (x)  the  first  (top)  and  second  (bottom) 
moments  of  the  probability  of  finding  hot  products  from 
the  lower  nozzle  with  an  isothermal  stream  in  the  upper 
nozzle. 


Fig.  9.  Examples  of  PIV  images  for  the  different  cases 
with  c/>  the  equivalence  ratio  of  the  JP-10  mixtures  in  the 
upper  stream.  Lines  indicate  detected  reaction  zones  and 
the  arrow  an  example  of  pocket  formation.  The  images 
are  inverted  for  clarity. 


mixing  layer  thickness  (<5,),  defined  as  the  region 
between  the  5th  and  95th  percentile,  was  esti¬ 
mated  from  c  and  the  ratio  bjLj  remains  effec¬ 
tively  unchanged  with  a  small  increase  from  3.2 
to  3.5  as  <p  is  reduced  from  0.6  to  0.4. 

Individual  images  with  reaction  zones  detected 
via  the  density  segregation  technique  were  further 
analysed  in  order  to  investigate  combustion  mode 
transitions.  The  mean  flame  lengths  in  a  radial 
interrogation  region  of  ±D/2  from  the  centreline 
were  calculated  and  an  integration  performed, 
assuming  axial  symmetry,  on  detected  reaction 
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Equivalence  Ratio  [-] 

Fig.  10.  Statistical  analyses  of  reaction  zones  from  PIV 
images  with  the  interface  area  obtained  by  assuming 
axial  symmetry  from  —  £>/ 2  to  D/2.  For  0  =  0  the 
stagnation  surfaces  were  detected  and  the  value  used  to 
normalise  results. 


zones  within  the  width  of  the  nozzle.  The  norma¬ 
lised  results  are  shown  in  Fig.  10  and,  as  suggested 
above,  a  distinct  transition  in  the  interface  (sur¬ 
face)  area  was  found  for  tj>  ^  0.4.  Turbulent  burn¬ 
ing  velocities  ( ST )  were  obtained  in  the  manner 
described  by  Lawn  and  Schefer  [28]  using  a  mean 
progress  variable  (c)  reference  point  of  0.02.  The 
resulting  values  are  shown  in  Table  1.  The  trend 
for  ST  follows  the  resolved  surface  area.  The  evo¬ 
lution  of  the  5/95  interface  thickness  with 

decreasing  equivalence  ratio  was  also  estimated 
from  the  PIV  images,  as  outlined  above,  with  val¬ 
ues  listed  in  Table  1. 

The  laminar  flame  properties  required  for  com¬ 
parisons  with  S1,  and  the  evaluation  of  the  corre¬ 
sponding  Da  numbers  are  not  readily  available 


Table  1 

Turbulent  burning  velocities  (Sf)  of  JP-10  and  the  local 
axial  rms  fluctuations  («')  at  locations  where  c  =  0.02 
[28]  along  the  centre  axis  of  the  burner.  Damkohler 
(Da  =  (Lr/Sf)/(urms/SL))  numbers  [16]  based  on  mea¬ 
sured  values  of  urms  =  0.69  m/s  and  L/=3.1mm.  The 
laminar  flame  parameters  (Sf  and  SL )  were  obtained 
from  computations  of  the  corresponding  cyclopentadi- 

ene  flames  [33].  The  estimated  interface  thickness  0?/)  is 
also  shown. 


<t> 

il 

[m/s] 

St 

[m/s] 

4 

[mm] 

5f 

[mm] 

SL 

[m/s] 

Da 

H 

0.0 

1.13 

- 

1.15 

- 

- 

- 

0.2 

1.18 

2.43 

1.57 

1.03 

0.05 

0.23 

0.4 

1.20 

2.56 

1.15 

0.94 

0.06 

0.28 

0.6 

1.14 

2.87 

0.78 

0.76 

0.14 

0.84 

0.8 

1.11 

3.34 

0.58 

0.34 

0.47 

6.0 

for  JP-10.  The  relationship  of  Parsinejad  et  al. 

[29]  for  JP-10  was  applied  at  400  K  to  estimate 
the  laminar  burning  velocity  ( SL )  at  0  =  0.80 
resulting  in  a  value  of  0.32  m/s.  By  contrast,  lam¬ 
inar  flame  calculations  performed  by  Courty  et  al. 

[30]  at  398  K,  using  the  JP-10  chemistry  of  Li  et  al. 

[31] ,  suggest  SL  =  0.51  m/s  at  the  same  stoichiom¬ 
etry.  Furthermore,  the  laminar  flame  thickness 
(Sf)  needs  to  be  estimated.  Davidson  et  al.  [32] 
suggested  that  cyclopentene  (C5J/8)  and  cyclo- 
pentadiene  ( C5H6 )  are  the  key  breakdown  prod¬ 
ucts  of  JP-10.  Accordingly,  laminar  flame 
calculations  featuring  the  more  reactive  C5H6 
were  performed  using  the  chemistry  of  Robinson 
and  Lindstedt  [33]  in  order  to  estimate  flame 
parameters  at  400  K.  The  computations  featured 
the  experimental  geometry  and  boundary  condi¬ 
tions  and  were  well  resolved  with  400  locally 
refined  control  volumes  using  the  standard  stag¬ 
nation  flow  co-ordinate  transformation  [34]  mod¬ 
ified  to  obtain  the  rate  of  strain  (a  ~  275/s)  as  an 
eigenvalue.  The  value  of  5f  was  determined  from 
the  computed  5/95  thickness  of  the  reaction  layer 
based  on  the  heat  release.  The  computed  heat 
release  and  temperature  based  reaction  progress 
variable  profiles  are  shown  in  Fig.  4.  The  thick¬ 
ness  of  the  reaction  layer  remains  comparatively 
close  to  the  experiment  despite  the  significant 
broadening.  The  corresponding  turbulent  back- 
to-back  opposed  jet  flames  extinguish  at  0  ~  0.7 
predominantly  due  to  the  impact  of  turbulent 
strain.  Hence,  while  a  conventional  turbulent 
flame  cannot  be  sustained,  combustion  remains 
supported  by  the  hot  combustion  products  in  a 
manner  that  is  consistent  with  the  HCDI  regime 


Fig.  11.  Axial  conditional  velocities  along  the  burner 
axis.  (O)  reactant  ( ur )  and  (□)  product  (up)  velocities  for 
JP-10  flames.  Top  left  <j>  =  0.20,  top  right  0  =  0.40, 
bottom  left  cf>  =  0.60  and  bottom  right  0  =  0.80.  Reac¬ 
tant  (ur)  velocities  with  an  isothermal  mixture  ( 0  =  0)  in 
the  upper  nozzle  are  also  shown  (x)  to  highlight 
the  increasing  deviation  with  an  increase  in  the 
stoichiometry. 
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Fig.  12.  Mean  axial  slip  velocities  us  (left)  and  scalar 
fluxes  u'c'  (right)  along  the  burner  axis.  Equivalence 
ratios  of  0.20  (O),  0.40  (O),  0.60  (A)  and  0.80  (□)  for 
lean  premixed  JP-10  flames. 

of  flameless  oxidation  [12].  The  relationship  to 
turbulent  combustion  regimes  [16,35]  can  be 
explored  further.  Velocity  fluctuations 
(urms  =  sj (ua  +  2t/2)/3  ~  0.69  m/s)  in  the  reac¬ 
tants  and  the  estimated  Damkohler  numbers  are 
shown  in  Table  1.  The  transition  from  the  corru¬ 
gated  flamelet  regime  to  distributed  reaction  zones 
is  estimated  to  occur  at  <j>  ~  0.8  and  the  further 
transition  to  the  well-stirred  regime  at  <j>  ~  0.6. 
The  density  segregation  technique  permits  a  de¬ 
tailed  analysis  of  conditional  statistics  of  interest 
to  the  evaluation  of  modelling  assumptions  such 
as  those  present  in  moment  [36]  and  transported 
PDF  [37]  based  approaches.  Mean  axial  reactant 
( ur )  and  product  ( up )  velocities  are  shown  in 
Fig.  1 1  and  conditional  velocities  derived  from 
density  segregation  can  be  used  to  determine  slip 
velocities  us,  as  defined  by  Bray  [38],  and  mean 
scalar  fluxes  u'c'  as  shown  in  Fig.  12.  Conditional 
velocities  and  scalar  fluxes  show  significant  devia¬ 
tions  from  the  isothermal  case  when  the  equiva¬ 
lence  ratio  was  raised  to  0.60  and  0.80.  Notably, 
scalar  fluxes  show  a  transition  to  counter-gradient 
transport  as  combustion  moves  towards  the  cor¬ 
rugated  flamelet  regime.  The  data  confirm  the 
continuous  transition  from  a  propagating  pre¬ 
mixed  turbulent  flame  to  flameless  combustion 
as  the  Da  number  is  reduced  with  the  overall 
structure  of  the  mixing  layer  increasingly  domi¬ 
nated  by  gradient  transport. 


5.  Conclusions 

The  continuous  transition  of  lean  premixed  tur¬ 
bulent  JP-10  flames  from  the  corrugated  flamelet 
regime  of  combustion  to  a  HCDI  related  regime 
of  flameless  combustion  has  been  investigated 
using  an  opposed  jet  geometry  with  hot  combus¬ 
tion  products  emerging  from  the  lower  nozzle. 
The  upper  nozzle  featured  JP-10/air  mixtures  with 
equivalence  ratios  of  0.20,  0.40,  0.60  and  0.80  as 
well  as  an  isothermal  (non-reacting)  case  (0  =  0) 
used  for  comparison  purposes.  The  observed  com¬ 
bustion  regime  transitions  were  quantified  using 
PIV  combined  with  a  density  segregation  tech¬ 
nique  for  the  determination  of  conditional  statis¬ 


tics.  The  growth  in  the  turbulent  burning  velocity 
and  flame  surface  area  was  also  quantified  for  the 
more  reactive  mixtures.  In  addition  to  the  move 
to  a  hot  combustion  product  supported  regime  of 
flameless  combustion,  it  was  also  noted  from  mea¬ 
sured  scalar  fluxes  that  turbulent  transport  changes 
from  a  gradient  to  a  counter-gradient  mode  when 
moving  into  the  corrugated  flamelet  regime. 
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